3.166 \(\int (e+f x)^2 \sin (a+b (c+d x)^2) \, dx\)

Optimal. Leaf size=256 \[ \frac{\sqrt{\frac{\pi }{2}} f^2 \cos (a) \text{FresnelC}\left (\sqrt{\frac{2}{\pi }} \sqrt{b} (c+d x)\right )}{2 b^{3/2} d^3}-\frac{\sqrt{\frac{\pi }{2}} f^2 \sin (a) S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right )}{2 b^{3/2} d^3}+\frac{\sqrt{\frac{\pi }{2}} \sin (a) (d e-c f)^2 \text{FresnelC}\left (\sqrt{\frac{2}{\pi }} \sqrt{b} (c+d x)\right )}{\sqrt{b} d^3}+\frac{\sqrt{\frac{\pi }{2}} \cos (a) (d e-c f)^2 S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right )}{\sqrt{b} d^3}-\frac{f (d e-c f) \cos \left (a+b (c+d x)^2\right )}{b d^3}-\frac{f^2 (c+d x) \cos \left (a+b (c+d x)^2\right )}{2 b d^3} \]

[Out]

-((f*(d*e - c*f)*Cos[a + b*(c + d*x)^2])/(b*d^3)) - (f^2*(c + d*x)*Cos[a + b*(c + d*x)^2])/(2*b*d^3) + (f^2*Sq
rt[Pi/2]*Cos[a]*FresnelC[Sqrt[b]*Sqrt[2/Pi]*(c + d*x)])/(2*b^(3/2)*d^3) + ((d*e - c*f)^2*Sqrt[Pi/2]*Cos[a]*Fre
snelS[Sqrt[b]*Sqrt[2/Pi]*(c + d*x)])/(Sqrt[b]*d^3) + ((d*e - c*f)^2*Sqrt[Pi/2]*FresnelC[Sqrt[b]*Sqrt[2/Pi]*(c
+ d*x)]*Sin[a])/(Sqrt[b]*d^3) - (f^2*Sqrt[Pi/2]*FresnelS[Sqrt[b]*Sqrt[2/Pi]*(c + d*x)]*Sin[a])/(2*b^(3/2)*d^3)

________________________________________________________________________________________

Rubi [A]  time = 0.339827, antiderivative size = 256, normalized size of antiderivative = 1., number of steps used = 11, number of rules used = 8, integrand size = 20, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.4, Rules used = {3433, 3353, 3352, 3351, 3379, 2638, 3385, 3354} \[ \frac{\sqrt{\frac{\pi }{2}} f^2 \cos (a) \text{FresnelC}\left (\sqrt{\frac{2}{\pi }} \sqrt{b} (c+d x)\right )}{2 b^{3/2} d^3}-\frac{\sqrt{\frac{\pi }{2}} f^2 \sin (a) S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right )}{2 b^{3/2} d^3}+\frac{\sqrt{\frac{\pi }{2}} \sin (a) (d e-c f)^2 \text{FresnelC}\left (\sqrt{\frac{2}{\pi }} \sqrt{b} (c+d x)\right )}{\sqrt{b} d^3}+\frac{\sqrt{\frac{\pi }{2}} \cos (a) (d e-c f)^2 S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right )}{\sqrt{b} d^3}-\frac{f (d e-c f) \cos \left (a+b (c+d x)^2\right )}{b d^3}-\frac{f^2 (c+d x) \cos \left (a+b (c+d x)^2\right )}{2 b d^3} \]

Antiderivative was successfully verified.

[In]

Int[(e + f*x)^2*Sin[a + b*(c + d*x)^2],x]

[Out]

-((f*(d*e - c*f)*Cos[a + b*(c + d*x)^2])/(b*d^3)) - (f^2*(c + d*x)*Cos[a + b*(c + d*x)^2])/(2*b*d^3) + (f^2*Sq
rt[Pi/2]*Cos[a]*FresnelC[Sqrt[b]*Sqrt[2/Pi]*(c + d*x)])/(2*b^(3/2)*d^3) + ((d*e - c*f)^2*Sqrt[Pi/2]*Cos[a]*Fre
snelS[Sqrt[b]*Sqrt[2/Pi]*(c + d*x)])/(Sqrt[b]*d^3) + ((d*e - c*f)^2*Sqrt[Pi/2]*FresnelC[Sqrt[b]*Sqrt[2/Pi]*(c
+ d*x)]*Sin[a])/(Sqrt[b]*d^3) - (f^2*Sqrt[Pi/2]*FresnelS[Sqrt[b]*Sqrt[2/Pi]*(c + d*x)]*Sin[a])/(2*b^(3/2)*d^3)

Rule 3433

Int[((g_.) + (h_.)*(x_))^(m_.)*((a_.) + (b_.)*Sin[(c_.) + (d_.)*((e_.) + (f_.)*(x_))^(n_)])^(p_.), x_Symbol] :
> Module[{k = If[FractionQ[n], Denominator[n], 1]}, Dist[k/f^(m + 1), Subst[Int[ExpandIntegrand[(a + b*Sin[c +
d*x^(k*n)])^p, x^(k - 1)*(f*g - e*h + h*x^k)^m, x], x], x, (e + f*x)^(1/k)], x]] /; FreeQ[{a, b, c, d, e, f,
g, h}, x] && IGtQ[p, 0] && IGtQ[m, 0]

Rule 3353

Int[Sin[(c_) + (d_.)*((e_.) + (f_.)*(x_))^2], x_Symbol] :> Dist[Sin[c], Int[Cos[d*(e + f*x)^2], x], x] + Dist[
Cos[c], Int[Sin[d*(e + f*x)^2], x], x] /; FreeQ[{c, d, e, f}, x]

Rule 3352

Int[Cos[(d_.)*((e_.) + (f_.)*(x_))^2], x_Symbol] :> Simp[(Sqrt[Pi/2]*FresnelC[Sqrt[2/Pi]*Rt[d, 2]*(e + f*x)])/
(f*Rt[d, 2]), x] /; FreeQ[{d, e, f}, x]

Rule 3351

Int[Sin[(d_.)*((e_.) + (f_.)*(x_))^2], x_Symbol] :> Simp[(Sqrt[Pi/2]*FresnelS[Sqrt[2/Pi]*Rt[d, 2]*(e + f*x)])/
(f*Rt[d, 2]), x] /; FreeQ[{d, e, f}, x]

Rule 3379

Int[(x_)^(m_.)*((a_.) + (b_.)*Sin[(c_.) + (d_.)*(x_)^(n_)])^(p_.), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplif
y[(m + 1)/n] - 1)*(a + b*Sin[c + d*x])^p, x], x, x^n], x] /; FreeQ[{a, b, c, d, m, n, p}, x] && IntegerQ[Simpl
ify[(m + 1)/n]] && (EqQ[p, 1] || EqQ[m, n - 1] || (IntegerQ[p] && GtQ[Simplify[(m + 1)/n], 0]))

Rule 2638

Int[sin[(c_.) + (d_.)*(x_)], x_Symbol] :> -Simp[Cos[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 3385

Int[((e_.)*(x_))^(m_.)*Sin[(c_.) + (d_.)*(x_)^(n_)], x_Symbol] :> -Simp[(e^(n - 1)*(e*x)^(m - n + 1)*Cos[c + d
*x^n])/(d*n), x] + Dist[(e^n*(m - n + 1))/(d*n), Int[(e*x)^(m - n)*Cos[c + d*x^n], x], x] /; FreeQ[{c, d, e},
x] && IGtQ[n, 0] && LtQ[n, m + 1]

Rule 3354

Int[Cos[(c_) + (d_.)*((e_.) + (f_.)*(x_))^2], x_Symbol] :> Dist[Cos[c], Int[Cos[d*(e + f*x)^2], x], x] - Dist[
Sin[c], Int[Sin[d*(e + f*x)^2], x], x] /; FreeQ[{c, d, e, f}, x]

Rubi steps

\begin{align*} \int (e+f x)^2 \sin \left (a+b (c+d x)^2\right ) \, dx &=\frac{\operatorname{Subst}\left (\int \left (d^2 e^2 \left (1+\frac{c f (-2 d e+c f)}{d^2 e^2}\right ) \sin \left (a+b x^2\right )+2 d e f \left (1-\frac{c f}{d e}\right ) x \sin \left (a+b x^2\right )+f^2 x^2 \sin \left (a+b x^2\right )\right ) \, dx,x,c+d x\right )}{d^3}\\ &=\frac{f^2 \operatorname{Subst}\left (\int x^2 \sin \left (a+b x^2\right ) \, dx,x,c+d x\right )}{d^3}+\frac{(2 f (d e-c f)) \operatorname{Subst}\left (\int x \sin \left (a+b x^2\right ) \, dx,x,c+d x\right )}{d^3}+\frac{(d e-c f)^2 \operatorname{Subst}\left (\int \sin \left (a+b x^2\right ) \, dx,x,c+d x\right )}{d^3}\\ &=-\frac{f^2 (c+d x) \cos \left (a+b (c+d x)^2\right )}{2 b d^3}+\frac{f^2 \operatorname{Subst}\left (\int \cos \left (a+b x^2\right ) \, dx,x,c+d x\right )}{2 b d^3}+\frac{(f (d e-c f)) \operatorname{Subst}\left (\int \sin (a+b x) \, dx,x,(c+d x)^2\right )}{d^3}+\frac{\left ((d e-c f)^2 \cos (a)\right ) \operatorname{Subst}\left (\int \sin \left (b x^2\right ) \, dx,x,c+d x\right )}{d^3}+\frac{\left ((d e-c f)^2 \sin (a)\right ) \operatorname{Subst}\left (\int \cos \left (b x^2\right ) \, dx,x,c+d x\right )}{d^3}\\ &=-\frac{f (d e-c f) \cos \left (a+b (c+d x)^2\right )}{b d^3}-\frac{f^2 (c+d x) \cos \left (a+b (c+d x)^2\right )}{2 b d^3}+\frac{(d e-c f)^2 \sqrt{\frac{\pi }{2}} \cos (a) S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right )}{\sqrt{b} d^3}+\frac{(d e-c f)^2 \sqrt{\frac{\pi }{2}} C\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right ) \sin (a)}{\sqrt{b} d^3}+\frac{\left (f^2 \cos (a)\right ) \operatorname{Subst}\left (\int \cos \left (b x^2\right ) \, dx,x,c+d x\right )}{2 b d^3}-\frac{\left (f^2 \sin (a)\right ) \operatorname{Subst}\left (\int \sin \left (b x^2\right ) \, dx,x,c+d x\right )}{2 b d^3}\\ &=-\frac{f (d e-c f) \cos \left (a+b (c+d x)^2\right )}{b d^3}-\frac{f^2 (c+d x) \cos \left (a+b (c+d x)^2\right )}{2 b d^3}+\frac{f^2 \sqrt{\frac{\pi }{2}} \cos (a) C\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right )}{2 b^{3/2} d^3}+\frac{(d e-c f)^2 \sqrt{\frac{\pi }{2}} \cos (a) S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right )}{\sqrt{b} d^3}+\frac{(d e-c f)^2 \sqrt{\frac{\pi }{2}} C\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right ) \sin (a)}{\sqrt{b} d^3}-\frac{f^2 \sqrt{\frac{\pi }{2}} S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right ) \sin (a)}{2 b^{3/2} d^3}\\ \end{align*}

Mathematica [A]  time = 1.82488, size = 151, normalized size = 0.59 \[ \frac{2 \sqrt{2 \pi } \text{FresnelC}\left (\sqrt{\frac{2}{\pi }} \sqrt{b} (c+d x)\right ) \left (2 b \sin (a) (d e-c f)^2+f^2 \cos (a)\right )+2 \sqrt{2 \pi } S\left (\sqrt{b} \sqrt{\frac{2}{\pi }} (c+d x)\right ) \left (2 b \cos (a) (d e-c f)^2-f^2 \sin (a)\right )-4 \sqrt{b} f (-c f+2 d e+d f x) \cos \left (a+b (c+d x)^2\right )}{8 b^{3/2} d^3} \]

Antiderivative was successfully verified.

[In]

Integrate[(e + f*x)^2*Sin[a + b*(c + d*x)^2],x]

[Out]

(-4*Sqrt[b]*f*(2*d*e - c*f + d*f*x)*Cos[a + b*(c + d*x)^2] + 2*Sqrt[2*Pi]*FresnelS[Sqrt[b]*Sqrt[2/Pi]*(c + d*x
)]*(2*b*(d*e - c*f)^2*Cos[a] - f^2*Sin[a]) + 2*Sqrt[2*Pi]*FresnelC[Sqrt[b]*Sqrt[2/Pi]*(c + d*x)]*(f^2*Cos[a] +
2*b*(d*e - c*f)^2*Sin[a]))/(8*b^(3/2)*d^3)

________________________________________________________________________________________

Maple [B]  time = 0.008, size = 669, normalized size = 2.6 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((f*x+e)^2*sin(a+(d*x+c)^2*b),x)

[Out]

-1/2*f^2/d^2/b*x*cos(b*d^2*x^2+2*b*c*d*x+b*c^2+a)-f^2*c/d*(-1/2/d^2/b*cos(b*d^2*x^2+2*b*c*d*x+b*c^2+a)-1/2*c/d
*2^(1/2)*Pi^(1/2)/(d^2*b)^(1/2)*(cos((b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelS(2^(1/2)/Pi^(1/2)/(d^2*b)^(1
/2)*(b*d^2*x+b*c*d))-sin((b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelC(2^(1/2)/Pi^(1/2)/(d^2*b)^(1/2)*(b*d^2*x
+b*c*d))))+1/4*f^2/d^2/b*2^(1/2)*Pi^(1/2)/(d^2*b)^(1/2)*(cos((b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelC(2^(
1/2)/Pi^(1/2)/(d^2*b)^(1/2)*(b*d^2*x+b*c*d))+sin((b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelS(2^(1/2)/Pi^(1/2
)/(d^2*b)^(1/2)*(b*d^2*x+b*c*d)))-e*f/d^2/b*cos(b*d^2*x^2+2*b*c*d*x+b*c^2+a)-e*f*c/d*2^(1/2)*Pi^(1/2)/(d^2*b)^
(1/2)*(cos((b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelS(2^(1/2)/Pi^(1/2)/(d^2*b)^(1/2)*(b*d^2*x+b*c*d))-sin((
b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelC(2^(1/2)/Pi^(1/2)/(d^2*b)^(1/2)*(b*d^2*x+b*c*d)))+1/2*2^(1/2)*Pi^(
1/2)/(d^2*b)^(1/2)*e^2*(cos((b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelS(2^(1/2)/Pi^(1/2)/(d^2*b)^(1/2)*(b*d^
2*x+b*c*d))-sin((b^2*c^2*d^2-d^2*b*(b*c^2+a))/d^2/b)*FresnelC(2^(1/2)/Pi^(1/2)/(d^2*b)^(1/2)*(b*d^2*x+b*c*d)))

________________________________________________________________________________________

Maxima [C]  time = 3.63412, size = 3950, normalized size = 15.43 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x+e)^2*sin(a+b*(d*x+c)^2),x, algorithm="maxima")

[Out]

-1/8*sqrt(pi)*(((-I*cos(1/4*pi + 1/2*arctan2(0, b)) - I*cos(-1/4*pi + 1/2*arctan2(0, b)) - sin(1/4*pi + 1/2*ar
ctan2(0, b)) + sin(-1/4*pi + 1/2*arctan2(0, b)))*cos(a) - (cos(1/4*pi + 1/2*arctan2(0, b)) + cos(-1/4*pi + 1/2
*arctan2(0, b)) - I*sin(1/4*pi + 1/2*arctan2(0, b)) + I*sin(-1/4*pi + 1/2*arctan2(0, b)))*sin(a))*erf((I*b*d*x
+ I*b*c)/sqrt(I*b)) + ((-I*cos(1/4*pi + 1/2*arctan2(0, b)) - I*cos(-1/4*pi + 1/2*arctan2(0, b)) + sin(1/4*pi
+ 1/2*arctan2(0, b)) - sin(-1/4*pi + 1/2*arctan2(0, b)))*cos(a) + (cos(1/4*pi + 1/2*arctan2(0, b)) + cos(-1/4*
pi + 1/2*arctan2(0, b)) + I*sin(1/4*pi + 1/2*arctan2(0, b)) - I*sin(-1/4*pi + 1/2*arctan2(0, b)))*sin(a))*erf(
(I*b*d*x + I*b*c)/sqrt(-I*b)))*e^2/(d*sqrt(abs(b))) - 1/2*(((e^(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + e^(-I*b
*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) + (I*e^(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) - I*e^(-I*b*d^2*x^2 - 2
*I*b*c*d*x - I*b*c^2))*sin(a))*abs(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2) + ((((2*I*sqrt(pi)*(erf(sqrt(I*b*d^2*x^2
+ 2*I*b*c*d*x + I*b*c^2)) - 1) - 2*I*sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) +
2*(sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) + sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*
d*x - I*b*c^2)) - 1))*sin(a))*b*c*d*x + ((2*I*sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) -
2*I*sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) + 2*(sqrt(pi)*(erf(sqrt(I*b*d^2*x^2
+ 2*I*b*c*d*x + I*b*c^2)) - 1) + sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*sin(a))*b*c^
2)*cos(1/2*arctan2(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2, 0)) + ((2*(sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x
+ I*b*c^2)) - 1) + sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) + (-2*I*sqrt(pi)*(er
f(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) + 2*I*sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^
2)) - 1))*sin(a))*b*c*d*x + (2*(sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) + sqrt(pi)*(erf(
sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) + (-2*I*sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x
+ I*b*c^2)) - 1) + 2*I*sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*sin(a))*b*c^2)*sin(1/2*
arctan2(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2, 0)))*sqrt(abs(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2)))*e*f/(b*d^2*abs(4
*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2)) + 1/32*(16*(b*d^2*x^2 + 2*b*c*d*x + b*c^2)^2*((e^(I*b*d^2*x^2 + 2*I*b*c*d*x
+ I*b*c^2) + e^(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) + (I*e^(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) -
I*e^(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*sin(a))*c + ((((I*sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I
*b*c^2)) - 1) - I*sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) + (sqrt(pi)*(erf(sqrt
(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) + sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*
sin(a))*b*c^2*d*x + ((I*sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) - I*sqrt(pi)*(erf(sqrt(-
I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) + (sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2))
- 1) + sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*sin(a))*b*c^3)*abs(4*b*d^2*x^2 + 8*b*c
*d*x + 4*b*c^2)*cos(1/2*arctan2(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2, 0)) + (((sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2
*I*b*c*d*x + I*b*c^2)) - 1) + sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) + (-I*sqr
t(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) + I*sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x -
I*b*c^2)) - 1))*sin(a))*b*c^2*d*x + ((sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2)) - 1) + sqrt(pi
)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*cos(a) + (-I*sqrt(pi)*(erf(sqrt(I*b*d^2*x^2 + 2*I*b*c
*d*x + I*b*c^2)) - 1) + I*sqrt(pi)*(erf(sqrt(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2)) - 1))*sin(a))*b*c^3)*abs(4
*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2)*sin(1/2*arctan2(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2, 0)) + (((-4*I*gamma(3/2,
I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + 4*I*gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) - 4*(gamma
(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*sin(a))*b*d^3*x
^3 + ((-12*I*gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + 12*I*gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*
b*c^2))*cos(a) - 12*(gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x -
I*b*c^2))*sin(a))*b*c*d^2*x^2 + ((-12*I*gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + 12*I*gamma(3/2, -I*
b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) - 12*(gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + gamma(3/2,
-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*sin(a))*b*c^2*d*x + ((-4*I*gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c
^2) + 4*I*gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) - 4*(gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x +
I*b*c^2) + gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*sin(a))*b*c^3)*cos(3/2*arctan2(4*b*d^2*x^2 + 8*b
*c*d*x + 4*b*c^2, 0)) - ((4*(gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + gamma(3/2, -I*b*d^2*x^2 - 2*I*b
*c*d*x - I*b*c^2))*cos(a) - (4*I*gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) - 4*I*gamma(3/2, -I*b*d^2*x^2
- 2*I*b*c*d*x - I*b*c^2))*sin(a))*b*d^3*x^3 + (12*(gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + gamma(3/
2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) - (12*I*gamma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) - 12*
I*gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*sin(a))*b*c*d^2*x^2 + (12*(gamma(3/2, I*b*d^2*x^2 + 2*I*b*
c*d*x + I*b*c^2) + gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) - (12*I*gamma(3/2, I*b*d^2*x^2 + 2
*I*b*c*d*x + I*b*c^2) - 12*I*gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*sin(a))*b*c^2*d*x + (4*(gamma(3
/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) + gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*cos(a) - (4*I*gam
ma(3/2, I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2) - 4*I*gamma(3/2, -I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2))*sin(a))*b
*c^3)*sin(3/2*arctan2(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2, 0)))*sqrt(abs(4*b*d^2*x^2 + 8*b*c*d*x + 4*b*c^2)))*f^
2/((b*d^2*x^2 + 2*b*c*d*x + b*c^2)^2*b*d^3)

________________________________________________________________________________________

Fricas [A]  time = 1.83911, size = 513, normalized size = 2. \begin{align*} \frac{\sqrt{2}{\left (\pi f^{2} \cos \left (a\right ) + 2 \, \pi{\left (b d^{2} e^{2} - 2 \, b c d e f + b c^{2} f^{2}\right )} \sin \left (a\right )\right )} \sqrt{\frac{b d^{2}}{\pi }} \operatorname{C}\left (\frac{\sqrt{2} \sqrt{\frac{b d^{2}}{\pi }}{\left (d x + c\right )}}{d}\right ) - \sqrt{2}{\left (\pi f^{2} \sin \left (a\right ) - 2 \, \pi{\left (b d^{2} e^{2} - 2 \, b c d e f + b c^{2} f^{2}\right )} \cos \left (a\right )\right )} \sqrt{\frac{b d^{2}}{\pi }} \operatorname{S}\left (\frac{\sqrt{2} \sqrt{\frac{b d^{2}}{\pi }}{\left (d x + c\right )}}{d}\right ) - 2 \,{\left (b d^{2} f^{2} x + 2 \, b d^{2} e f - b c d f^{2}\right )} \cos \left (b d^{2} x^{2} + 2 \, b c d x + b c^{2} + a\right )}{4 \, b^{2} d^{4}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x+e)^2*sin(a+b*(d*x+c)^2),x, algorithm="fricas")

[Out]

1/4*(sqrt(2)*(pi*f^2*cos(a) + 2*pi*(b*d^2*e^2 - 2*b*c*d*e*f + b*c^2*f^2)*sin(a))*sqrt(b*d^2/pi)*fresnel_cos(sq
rt(2)*sqrt(b*d^2/pi)*(d*x + c)/d) - sqrt(2)*(pi*f^2*sin(a) - 2*pi*(b*d^2*e^2 - 2*b*c*d*e*f + b*c^2*f^2)*cos(a)
)*sqrt(b*d^2/pi)*fresnel_sin(sqrt(2)*sqrt(b*d^2/pi)*(d*x + c)/d) - 2*(b*d^2*f^2*x + 2*b*d^2*e*f - b*c*d*f^2)*c
os(b*d^2*x^2 + 2*b*c*d*x + b*c^2 + a))/(b^2*d^4)

________________________________________________________________________________________

Sympy [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \left (e + f x\right )^{2} \sin{\left (a + b c^{2} + 2 b c d x + b d^{2} x^{2} \right )}\, dx \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x+e)**2*sin(a+b*(d*x+c)**2),x)

[Out]

Integral((e + f*x)**2*sin(a + b*c**2 + 2*b*c*d*x + b*d**2*x**2), x)

________________________________________________________________________________________

Giac [C]  time = 1.23698, size = 952, normalized size = 3.72 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x+e)^2*sin(a+b*(d*x+c)^2),x, algorithm="giac")

[Out]

1/4*I*sqrt(2)*sqrt(pi)*erf(-1/2*sqrt(2)*sqrt(b*d^2)*(-I*b*d^2/sqrt(b^2*d^4) + 1)*(x + c/d))*e^(I*a + 2)/(sqrt(
b*d^2)*(-I*b*d^2/sqrt(b^2*d^4) + 1)) - 1/4*I*sqrt(2)*sqrt(pi)*erf(-1/2*sqrt(2)*sqrt(b*d^2)*(I*b*d^2/sqrt(b^2*d
^4) + 1)*(x + c/d))*e^(-I*a + 2)/(sqrt(b*d^2)*(I*b*d^2/sqrt(b^2*d^4) + 1)) - 1/2*(I*sqrt(2)*sqrt(pi)*c*f*erf(-
1/2*sqrt(2)*sqrt(b*d^2)*(-I*b*d^2/sqrt(b^2*d^4) + 1)*(x + c/d))*e^(I*a + 1)/(sqrt(b*d^2)*(-I*b*d^2/sqrt(b^2*d^
4) + 1)) + f*e^(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2 + I*a + 1)/(b*d))/d - 1/2*(-I*sqrt(2)*sqrt(pi)*c*f*erf(-1/
2*sqrt(2)*sqrt(b*d^2)*(I*b*d^2/sqrt(b^2*d^4) + 1)*(x + c/d))*e^(-I*a + 1)/(sqrt(b*d^2)*(I*b*d^2/sqrt(b^2*d^4)
+ 1)) + f*e^(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2 - I*a + 1)/(b*d))/d - 1/8*(-I*sqrt(2)*sqrt(pi)*(2*b*c^2*f^2
+ I*f^2)*erf(-1/2*sqrt(2)*sqrt(b*d^2)*(-I*b*d^2/sqrt(b^2*d^4) + 1)*(x + c/d))*e^(I*a)/(sqrt(b*d^2)*(-I*b*d^2/s
qrt(b^2*d^4) + 1)*b) + 2*I*(d*f^2*(-I*x - I*c/d) + 2*I*c*f^2)*e^(I*b*d^2*x^2 + 2*I*b*c*d*x + I*b*c^2 + I*a)/(b
*d))/d^2 - 1/8*(I*sqrt(2)*sqrt(pi)*(2*b*c^2*f^2 - I*f^2)*erf(-1/2*sqrt(2)*sqrt(b*d^2)*(I*b*d^2/sqrt(b^2*d^4) +
1)*(x + c/d))*e^(-I*a)/(sqrt(b*d^2)*(I*b*d^2/sqrt(b^2*d^4) + 1)*b) + 2*I*(d*f^2*(-I*x - I*c/d) + 2*I*c*f^2)*e
^(-I*b*d^2*x^2 - 2*I*b*c*d*x - I*b*c^2 - I*a)/(b*d))/d^2